Superconductivity in compressed iron: Role of spin fluctuations 
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The recent discovery of superconductivity in hexagonal iron under pressure poses a question 
about whether it is of conventional (phonon) or unconventional (magnetic?) origin. We present 
first-principles calculations of the electron-phonon coupling in iron at P > 15 GPa, and argue 
that a conventional mechanism can explain the appearance of superconductivity, but not its rapid 
disappearance at P > 30 GPa. We suggest that spin fluctuations, ferro- and/or antiferromagnetic, 
play a crucial role in superconductivity in this case. 
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The rjecent report of superconductivity in Fc under 
pressuretl was among the latest news in the chain of dis- 
coveries of interesting cases of supercwnductivity in con- 
ventional materials: MgB2H, MgCNial, ZiZnM, and oth- 
ers. At first glance, there is nothing unexpected in the 
fact that iron, a transition metal, becomes superconduct- 
ing in a nonmagnetic phase. However, a closer look at 
this case uncovers a number of strange facts which do not 
square with this simplistic scenario. First, in the phase 
diagram (Fig.l) according to Rcf. |l| the superconductiv- 
ity appears at the same pressure (P « 15 GPa), at which 
the bcc-hcp phase transition takes place. In other words, 
the bcc phase is not superconducting (and is magnetic), 
while the hep phase is already superconducting at the 
lowest pressure at which it is stable. What is surprising, 
however, is that the critical temperature is zero at the 
phase transition pressure, and then grows gradually with 
pressure. If the superconductivity were due to phonons, 
its appearence exactly at the crystallographic phase tran- 
sition would be just a coincidence, since there is no reason 
for superconductivity to be suppressed near such a (first 
order) phase transition. Second, the superconductivity 
disappears very rapidly with increasing pressure. Essen- 
tially, superconductivity occurs only in a narrow range of 
compression, 132 Bohr^< V < 145 Bohr"^, where V is the 
unit cell volume (with two atoms per cell). As we show 
below, neither electronic (e.g., density of states) char- 
acteristics nor phonon frequencies change with pressure 
rapidly enough to reduce the critical temperature from 
nearly 2 K to zero at such a small relative compression. 

On the other hand, there are cases in nature when su- 
perconductivity exists only in a narrow pressure range, as 
in the_recently discovered ferromagnetic superconductor 
UGe2.El The accepted explanation for such behavior in 
this case is that superconductivity there is believed to be 
induced by spin fluctuations, which, first, stiffen rapidly 
away from the quantum critical point (that is, the pres- 
sure that corresponds to the Curie temperature Tc — 0), 
and, second, grow in magnitude close to this point. Thus, 
the critical temperature goes to zero at the critical point, 
then rises away from it, but disappears again whftH the 
energy of the spin fluctuations becomes too highJj'El. 

Finally, it is not even clear whether hep iron is in- 
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deed nonmagnetic over its whole range of stability. Den- 
sity functional calculations using the Generalized Gra- 
dient Approximation (GGA), which is usually very re- 
liable for 3d-.metals, clearly show an antiferromagnetic 
ground stateH, or a noncoUinear magnetic stateH, but not 
a paramagnetic one. Furthermore, the calculated elastic 
properties in magnetic hep iron are in good agreement 
with experiment, while those for a. nonmagnetic iron se- 
riously disagree with experiments Admittedly, there is 
one experiment that is hard to explain if hep iron is mag- 
netic even, locally, namely Mossbauer spectroscopy under 
pressurecj. Given this contradiction, we take the point 
of view that the question of magnetism in hep iron is still 
open, and even if the hep iron is in fact non-magnetic, it 
should be very close to a magnetic instability. 

In 1979, well before superconductivity was discovered 
in hep iron. Wok 
may play a role 
sures at which the hep iron is stable, it must be very close 
to an antiferromagnetic instability, which should lead to 
pair-breaking spin fluctuations. He also pointed out that 
these fluctuations should disappear at higher pressure, as 
the density of states becomes smaller, and thus at some 
pressure superconductivity is to be expected. Of course, 
what is missing from this scenario is the fact that super- 
conductivity disappears again rather soon when pressure 
is increased. 

In order to elucidate the potential for conventional, 
electron-phonon superconductivity in iron we performed 
first-principles Linear Augmented Plane Wave (LAPW) 
calculations in GGA of the following quantities: (a) the 
electron-ion matrix element, W, (b) the density of states 
at the Fermi level, N, and (c) the Stoner factor, /. Su- 
perconductivity in transition metals has been subject to 
numerous studies in the 70's and 80's (see, e.g., Ref. p^for 
review), and is quantitatively understood. The McMillan 
formula. 
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describes reasonably well. The parameters of the for- 
mula have the following meanings: wiog is the logarith- 
mically averaged phonon frequency, /i* is the Coulomb 
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pseudopotential, which for transition metals is usually 
close to 0.1, and A is the electron-phonon coupling con- 
stant, A = r]M^^ {^^^) 1 where r/ = NW is the so-called 
Hopfield parameter, M is the ion mass, and is the 

average inverse square of the phonon frequency. Note 
that the electron-phonon coupling is responsible for both 
pairing interaction and mass renormalization. The con- 
stant A measures both effects and appears three times in 
the McMillan formula: once as a measure of the pairing 
strength, and twice (in 1 -I- A and in 1 -f 0.54A) because 
of mass renormalization. 

In at least one transition metal, Pd, the electron- 
phonon coupling is definitely strong enough to render 
the material superconducting but, because it is near a 
magnetic instability, superconductivity is suppressed. In 
other words, charge fluctuations (phonons) coexist in Pd 
with spin- fluctuations (paramagnons), which have a pair- 
breaking effect for s-wave symmetry, but contribute to 
the mass renormalization like phonons. These effects can 
be approximately accounted for by modifying the McMil- 
lan formula as follows: 
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where the subscripts stand for " phonons" and " spin fluc- 
tuations" , so that Xsf is an electron-paramagnon cou- 
pling constant. 

Since the Migdal theorem does not hold for param- 
agnons, it is hardly possible to derive a formula for A^/ 
and compute it from first principles. Within the frame- 
work of the Stoner model, and near the magnetic insta- 
bility, one can jestimate coupling with ferromagnetic spin 
fiuctuations aa3 
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where / is the Stoner factor and a characterizes the 
band structure andj-tJie (7— dependence of the bare mag- 
netic susceptibilitjxJ. The most accurate way to es- 
timate / is by performing fixed spin moment calcula- 
tions for small moments and taking the second deriva- 
tive of the total energy with respect to the fixed moment, 
/ = N^^ - 2(PE/d'W?. With this method we calculated 
/ =0.075 Ry/atom, using fixed-spin-moment LAPW to- 
tal energies aXV— 146.53 Bohr"^ . Note that / is practi- 
cally pressure-independent. 

For antiferromagnetic spin fiuctuations EqJ^ is not ap- 
plicable, but instead one can write, as a first approxima- 
tion. 
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where x is the noninteracting spin susceptibility at the 
wave vector corresponding to the antiferromagnetic or- 
dering. The main effect of pressure on the electronic 



structure of hep iron is an overall scaling of the band- 
width. Thus, one can assume that x ~ f^N, where /3 is a 
pressure independent constant. 

All quantities entering Eq. |2[ except a and /?, were cal- 
culated from first principles, using the approximations 
outlined above with the phonon frequencies, calculated 
elsewhcreo, which agree well with cxperimentll3. It is 
worth noting that because of all these approximations 
we do not expect to obtain quantitatively the value for 
Tc, but we do expect to describe correctly the trends as- 
sociated with the pressurell3. 

In the calculations of the electron-phonon Hopfield pa- 
rameter 77, we use|d-|the " rigid- muffin-tin" theory of Gas- 
pari and GyorffyJlJ According to this theory, which is 
known to work quite well for transition metals (see, e.g. 
Ref. ^,V = Y.i W^NiNi+i/N, where Wi are computed 
from particular integrals involving the radial wave func- 
tions, as defined in Ref. |l^, and Ni is the partial den- 
sity of states at the Fermi level with angular momentum 
I. We used self-consistent MT potentials with touching 
MT spjheres, calculated by the LAPW method. It is well 
knownll3 that for transition metals the main contribution 
to r] comes from I = 2, thus r] « WjNdNf/N « WjNf, 
as the ratio Nd/N is close to unity in transition metals. 
In this connection, it is worth noting that the volume de- 
pendence of 77 (cf. Fig.^ is totally different from that of 
the total density of states iV; in fact it has the opposite 
variation. 

The pressure-independent constants a and (3 were ad- 
justed so as to have the magnetic instability {(3NI — 1) 
near the hcp-bcc phase boundary, V « 145 Bohr^, and 
to have the maximal critical temperature of the right 
order. We used a — 0.13 and /? = 1.215. Qualitative de- 
pendence of Tc on pressure is not sensitive to the actual 
values of these constants. 

The results of our calculations are summarized in Table 
1, and compared with experiment on Fig.0. We observe 
that the phase diagram can be described qualitatively 
within the framework of the conventional scenario (pair- 
ing phonons plus pair-breaking spin fiuctuations), in the 
sense that the critical temperature rises with pressure 
near the bcc-fcc transition pressure, reaches a maximum, 
and then decays. The underlying physics is as follows: 
spin fluctuations, which get stronger with lowering pres- 
sure, suppress superconductivity, and at very high pres- 
sures the lattice stiffens and Tc goes down. However, 
the range of pressures at which superconductivity exists 
is much larger in the calculations. This is not acciden- 
tal; the only mechanism that prevents superconductivity 
at higher pressure is stiffening of the lattice, since the 
electronic part, 77, actually grows with pressure (Fig||). 
Indeed, hep Fe does stiffen rapidly under pressure. How- 
ever, it is by far not strong enough an effect to destabilize 
superconductivity at P > 35 GPa (corresponding to the 
atomic volume of 132 Bohr'^, which is only a 4% com- 
pression from the maximum Tc volume). In other words, 
elastic properties do not change fast enough with pres- 
sure to explain such an extraordinarily narrow range of 
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superconductivity. 

On the other hand, magnetic properties, at least as cal- 
culated within density functional theory, do change in the 
relevant range. Indeed, at V = 180 Bohr^, both antifer- 
romagnetic arrangements considered in Ref. ^ are stable 
and produce the sarae, magnetic moment, on the order 
of 2 /iB per Fe atom£3. At V = 145 Bohr^, the largest 
volume at which hep iron exists, both arrangements are 
still stable, and the corresponding magnetic moments are 
similar, although they reduce to values of « 0.5 ^b- One 
can interpret this as local magnetism on Fe with nearest- 
neighbor antiferromagnetic exchange (J). The exchange 
energy for the magnetic structure afml (we use the no- 
tation of Ref. ^) is zero, for afmll it is -2J per atom, 
and for a ferromagnetic state it is +6J per atom. Ap- 
parently, the latter energy is sufhciently large to render 
the ferromagnetic state unstable. In fact, and this has 
not been noticed before, while the ferromagnetic state 
is energetically unfavorable compared with the nonmag- 
netic state for all volumes up to = 157 Bohr'^, there is 
a metastable ferromagnetic state (metamagnetism) with 
M w 2.6 /Lts/atom, for all V > 145 Bohr^ (Fig.|) On the 
other hand, while antiferromagnetic states, stabilized by 
the exchange energy, are lower in energy than the non- 
magnetic one, the equilibrium magnetization is rather 
small, on the order of 1 fis- We find, in agreement with 
Ref. ^, that the antiferromagnetism disappears only at 
V « 120 Bohr3. 

While the picture derived from these calculations is at 
odds with that deduced from the Mossbauer experiment, 
in the sense that the calculated tendency to magnetism 
appears to be too strong, the main conclusion can be ex- 
pressed qualitatively as follows: the spectrum of the spin 
fluctuations changes with pressure rapidly and drastically 
near P « 20 GPa, and probably in a very nontrivial way, 
reflecting two different types of magnetic excitations: an- 
tiferromagnetic, and also ferromagnetic (the latter, pos- 
sibly, of highly nonlinear metamagnetic character). This 
strongly suggests an exclusive role played by spin fluc- 
tuations in iron superconductivity around P « 20 GPa. 
Whether the superconductivity is s-wave, induced by the 
electron-phonon interaction, with spin fluctuations be- 
ing ordinary pair breakers, or it is d-wave or p-wave, 
with spin fluctuations being the pairing agents, cannot 
be firmly established now, and should be clarified by fur- 
ther experimental studies. However, the latter scenario 
seems to be not only more exciting, but also more likely, 
if we believe that, according to the Mossbauer results, 
there are no localized magnetic moments in the whole su- 
perconductivity range, because this means that the spin 
fluctuations become weaker away from the bcc-hcp tran- 
sition, and their pair-breaking effect can only diminish 
with pressure. On the other hand, if superconductivity 
is d- or p-wave, and the spin fluctuations are pairing, the 
phase diagram looksjponsistent with the standard predic- 
tion for this scenarioQ. A key experiment in this situation 
would be measuring the effect of nonmagnetic impurities 
on superconductivity. 
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FIG. 1. Volume dependence of experimental and calcu- 
lated superconducting properties of hep iron. The elec- 
tron-phonon coupling constant. A, and the critical tem- 
perature for the s-wave superconductivity, Tc, (according 
to Eq.|2|) were computed as described in the text. The 
constant a, controlling the strength of the pair-breaking 
effect of spin fluctuations, was adjusted so as to have su- 
perconductivity disappear aX V k, 145 Bohr'', as in the 
experiment. The experimental critical temperatures from 
Ref. |l| are also plotted. For translating the pressure scale 
into tha-jj-olume scale we used the experimental equation 
of stateEj, which is also shown in the Figure. We also 
plotted the Stoner product, /TV. 




FIG. 2. Hopfield factor r] in Ry/Bohr'^ compared with 
the density of states at the Fermi level, A''. 
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FIG. 3. LAPW flxed spin moment calculations for fer- 
romagnetic hep iron. Magnetic moments and energies are 
given on a per cell (two atoms) basis. The curves corre- 
spond, from up down, to the volumes per atom: 169.97 
Bohr^ 160.44 Bohr^ 155.81 Bohr^ 151.27 Bohr^ and 
146.83 Bohr^ 



TABLE I. 
volumes. 



Calculated parameters at three representative 



V", Bohr-' 
P, GPa^"' 
?7, Ry/Bohr^ 



106.83 125.63 146.53 



2), cm" 

-1(6) 



-1(6) 



cjiog, cm 
-/V, st./Ry atom spin 
1/(1 - IN) 



K, Eq 
K, Eq 



2 (a = 0.013, /3 = 1.215) 



79 


48 


12 


0.185 


0.132 


0.099 


332 


260 


206 


343 


260 


200 


6.7 


8.3 


10.4 


2.0 


2.7 


4.6 


2.5 


4.1 


6.0 


1.0 


1.4 


0.0 



'^''^ These values are calculated from the data in Ref. |l^, 
assuming that the phonon density of states is uniformly 
distributed between the longitudinal and the transverse 
peaks. In fact, any reasonable assumption about the 
shape of the phonon DOS leads to very similar results. 
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